Ionization of clusters in strong X-ray laser pulses 
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The effect of intense X-ray laser interaction on argon clusters is studied theoretically with a mixed 
quantum/classical approach. In comparison to a single atom we find that ionization of the cluster 
is suppressed, which is in striking contrast to the observed behavior of rare-gas clusters in intense 
optical laser pulses. We have identified two effects responsible for this phenomenon: A high space 
charge of the cluster in combination with a small quiver amplitude and derealization of electrons 
in the cluster. We elucidate their impact for different field strengths and cluster sizes. 
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The advent of femtosecond laser pulses has triggered 
numerous activities in the field of laser-matter interac- 
tion. In particular novel, non-linear processes in atoms 
induced by the available high intense fields from these 
ultrashort pulses have been a challenge for experiment 
and theory [[l] . To understand more generally the energy 
transfer from intense laser light to matter, complex tar- 
gets like atomic clusters have been studied || |) . 
These experiments received considerable attention due 
to the observed dramatic effects like emission of very fast 
ions |^] and electrons || or the production of coherent 
X-ray radiation 0]. So far such studies are almost ex- 
clusively restricted to visible or infrared wavelengths || . 
The new X-ray free electron laser (xfel) sources, under 
construction at DESY in Hamburg |?J and at the LCLS 
in Stanford ||], will change this situation. They can de- 
liver intense laser pulses at high frequencies (from VUV 
to hard X-ray) and thus open a new regime of strong-field 
atomic physics. 

Here we present theoretical investigations of X-ray 
(frequency fiw = 350 eV, i.e. wavelength A f=a 3.5 nm) 
laser interaction with argon clusters at intensities of 
I « 3.5-10 14 ...3.5-10 18 W/cm 2 . In this laser regime 
the interaction is notably different from long-wavelength 
pulses which is evident from the ponderomotive energy 
£-pond ~ I/uj 2 . It represents the average kinetic energy 
of a free electron in a laser field, while Aa; ~ \/l /uj 2 is 
the quiver amplitude of the electron, i. e., the spatial ex- 
cursion in the laser field. For the given laser parameters 
one finds -Epond ~ 0.4 meV. . . 4eV, vastly different from 
-Epond ~ 20 eV . . . 200 keV for a 780-nm-laser at the same 
intensities. The small ponderomotive energy has two im- 
portant consequences: (i) Despite the high intensities the 
laser-atom interaction is of non-relativistic and pertur- 
bative nature. The latter is also clear from the so-called 
Keldysh parameter 7 which gives the ratio between the 
tunneling time and the laser period |^] and can be rewrit- 
ten in terms of the binding energy -Ebind and the pondero- 
motive energy E pond as 7 = (£'bind/2Epond) 1/2 - For the L 
and M-shells of argon with Sbind = 326, 249, 29.3, 15.8 eV 
we always find 7 > 1 and for the inner shells even 
7 ^> 1. Obviously the laser period is far too small 
for field- or tunnel-ionization. Rather, ionization is due 
to single-photon absorption, (ii) Ponderomotive effects 



are completely negligible with quiver amplitudes of 
Ax w 0.0003... 0.03 A. 

Under the aforementioned laser conditions one expects 
exhaustive ionization leading subsequently to a gigan- 
tic Coulomb explosion of the irradiated argon clusters, 
since ionization proceeds fundamentally different from 
that under optical laser impact. Ionization starts from 
the inside because photoionization cross sections at X-ray 
wavelengths are considerably higher for the inner shells 
than for the valence shells pj| . Typically the inverse 
rates are with 1 ... 10 fs |Q much smaller than the pulse 
length of about 100 fs. Hence, multiple single-photon ion- 
ization is possible, in particular because the inner-shell 
holes created by photoionization are refilled by Auger-like 
processes. The Auger decay is more or less independent 
of the atomic charge state and occurs fast, typical times 
are 0.2 . . . 5fs Jf|. Due to this almost instantaneous re- 
filling of the inner shells they can be ionized many times 
during the pulse and thus the atoms can be efficiently 
"pumped dry". This occurs "inside-out" and is the ex- 
act opposite to the ionization mechanism in the visible 
wavelengths regime where the electrons are removed from 
the outside like shells of an onion jl4j . Finally, local ionic 
charges in the cluster, generated in the course of the ion- 
ization process, may be screened by the weakly bound 
valence electrons of cluster atoms. Therefore, these elec- 
trons are allowed to tunnel through the barriers between 
the mother atom and neighboring ions. Such electrons 
will predominantly move inside the cluster since there 
is no field (such as the quasi-static electric field of low- 
frequency lasers) which would drive them away from the 
charged cluster into the continuum. 

In order to gain insight into the laser induced dynam- 
ics we use a mixed quantum/classical approach similar to 
those successfully applied in studying rare-gas clusters 
in strong pulses of visible light |lj, |l5[ |l(|. Electrons 
initially bound to atoms in the cluster are not explic- 
itly treated, they only enter through their binding en- 
ergy. The ionization from the mother atom in the clus- 
ter, also called inner ionization, is modelled statistically 
by sudden transitions according to quantum mechanical 
rates. An inner ionization event gives birth to a quasi- 
free electron which is subsequently propagated classically 
along with the ions and other quasi-free electrons with all 
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Coulomb forces included. Note that a quasi-free electron 
is not necessarily ionized with respect to the cluster as a 
whole, quasi-free means only inner-ionized with respect 
to the mother atom or ion. 

The ionization/excitation dynamics in strong X-ray 
pulses is considerably more complicated than field []1 4f| 
or tunnel jl^, [l6| ionization in optical pulses because 
inner-shell electrons and intra-atomic processes are in- 
volved. For the photoionization rates we use a parameter- 
ization of the cross section which covers all charge states 
of atomic argon jlTj . For high ionic charges the electronic 
binding energy can become so large that single photoion- 
ization is impossible. Non-dipole effects do not have any 
crucial impact apart from distortions of the angular dis- 
tribution of the photoelectrons [|l8) . For modelling the 
intra-atomic decay cascades we use direct Hartree-Fock 
calculations fUj . They provide branching ratios and de- 
cay rates for Auger, Coster-Kroning, and radiative tran- 
sitions. The latter ones are typically at least one order of 
magnitude slower Q . For the intra-cluster charge equi- 
libration we calculate tunnel (or over-the-barrier) rates 
in analogy to (l^, 111]. All the rates are used to decide 
at every time step of the simulation whether a particular 
transition occurs or not. Once the electrons are created 
with their respective kinetic energy (to guarantee energy 
conservation in the photo-absorption process or in the 
Auger decay jl9| ) they are propagated in the field of the 
other particles. 

From the physics described so far it is clear that the 
initial coupling to the laser is individual atom-photon in- 
teraction only. Hence, it is particularly interesting to see 
how the cluster environment affects the laser interaction 
which can be most clearly seen from calculations com- 
paring atoms and clusters. The laser pulse is given by 
f(t) = f ■ s(t) ■ sin(cji) with the electric field strength /, 
the laser frequency ui, and the pulse shape function 
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where T is the pulse length. Figure [l] shows the fi- 
nal charge per atom for two different clusters, Ari3 and 
Arss, compared to a single atom for a laser pulse with 
hu = 350 eV, T = 100 fs, and / = 0.1 . . . 10 au (atomic 
unit). At lower field strengths (/ = 0.1 au), only a single 
photoionization event with one subsequent Auger decay 
per atom occurs independently of the cluster size. Af- 
ter a quite steep rise the final charge starts to saturate 
in the atomic case for stronger fields (/ > 0.3 au) due 
to the fact that single photoionization becomes impossi- 
ble beyond a certain charge state of the ion at the given 
photon frequency. The rise of the final charge per atom 
in the cluster is notably weaker. Therefore — in strik- 
ing contrast to optical laser-cluster interaction — clusters 
are less effectively ionized at high fields than atoms. The 
reduction is even more pronounced for the larger cluster 
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FIG. 1: Average charge per atom for two cluster sizes Ari3 
and Ar 55 produced by an xfel pulse {Tiuj = 350 eV, T = 
100 fs) as a function of the field strength / compared to atomic 
argon as target. The points are averages from 10 simulations. 



(cf. Fig. |). One reason for the reduced ionization in the 
cluster compared to the atom is the much larger space 
charge produced in a cluster. At / = 10 au the total clus- 
ter charges Q after the pulse are Qi 7, Qi 3 w 65, and 
Q 55 ss 220 for atoms and clusters of 13 and 55 atoms, 
respectively. Such high space charges suppress ioniza- 
tion, because the absorption of one single photon trans- 
fers only a fixed amount of energy to the electron. In 
our model, these post-collisional interaction effects are 
taken into account by the propagation of the photoion- 
ized electrons. For larger cluster and higher fields, i. e. 
higher space charges, one finds that an increasing num- 
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FIG. 2: Absorbed energy per atom as a function of the field 
strength / for the three targets and the laser pulse shown in 
Fig. 
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ber of electrons are bound at the end of the pulse to one 
of the fragment ions. 

Figure |^ shows the energy absorbed from the laser 
pulse for the same clusters. This observable provides 
direct insight into the photoionization since all the other 
processes (intra-atomic decay, intra-cluster screening) are 
not influenced by the laser due to its high frequency. Sur- 
prisingly, for all field strengths considered, the absorption 
of photons itself is reduced in the environment provided 
by the cluster. This was unexpected because of the fact 
that predominantly deep-lying inner-shell electrons are 
affected and possible effects on the photoionization rates 
of these strongly localized electrons from neighboring ions 
are unlikely and in fact not contained in the model. 

The time evolution of the cluster dynamics reveals that 
a derealization of the valence electrons is indirectly re- 
sponsible. It has mainly two effects: Firstly, the pho- 
toionization cross sections become very small since the 
electrons are far away from the nucleus. Secondly, also 
the Auger decay rates are reduced because the overlap 
with the core holes becomes smaller. Before quantifying 
these effects we will look at the electron dynamics in de- 
tail. Figure |^ shows the time evolution of an Ar^ cluster 
for the same pulse parameters as before at a field strength 
of / = 1 au. At every time step each electron is assigned 
to the cluster ion with which it has the largest instant 
binding energy. In this way we define a "hopping" time 
which is the time an electron stays at the same ion before 
it moves to another one. The upper panel of Fig. || shows 
that the average "hopping" time is less than 1 fs during 
the first half of the pulse (t < 50 fs). This is very short 
compared to typical inverse ionization rates and indicates 
that the electrons move almost freely inside the cluster 
volume. Hence, the low barriers between the cluster ions 
in the early phase of the pulse make inner ionization very 
effective. As can be seen in the middle panel of Fig. ||, a 
considerable fraction of the quasi-free electrons has ener- 
gies above these barriers. This will change only towards 
the end of the pulse (t = 100 fs) when due to the expan- 
sion of the cluster (cf. the cluster radius in Fig. ||) the 
barriers rise again. 

It is difficult to assign photoionization or autoioniza- 
tion rates to this "sea" of electrons. However, one can 
give estimates of these rates using the fact that they de- 
pend in a characteristic manner on the energy of the re- 
spective bound electron. The quasi-free electrons have 
a binding energy Ei Tee (since they are still bound with 
respect to the full cluster). This energy will be com- 
pared to -©bound the binding energy of the weakest bound 
electron in the atom, which is among the bound elec- 
trons the most stable one against ionization. The ratios 
_Ef ree /i?bound, averaged over all the electrons, are shown 
in the lower panel of Fig. ||. Compared to the low rates of 
the bound electrons the photoionization rate falls off like 
(.Ef ree /i?bound) 7 ^ 2 (llj and the Auger decay rate roughly 
like (©f re e /©bound Due to this scaling the rates 

are reduced during the pulse by a factors of about 5 and 
20, respectively. Obviously, also those electrons with en- 



ergies below the barriers are still well above the highest 
bound electrons and hence fairly delocalized. One has to 
emphasize that this estimate is an upper bound since the 
rates decrease even more for higher angular momentum 
states O, 20 which are likely to be populated by the 
intra-cluster dynamics. Therefore, we regard it as save 
to neglect absorption of photons or autoionization of the 
quasi- free electrons. 

As we have seen in addition to simply lower absorption 
rates the reduction of photon absorption is also due to 
the fact that the inner shells to be ionized are no longer 
efficiently refilled by inter-atomic decay, i. e. the cluster is 
temporarily hollow. In order to assess the relative impor- 
tance of this effect compared to the suppression of ion- 




FIG. 3: Time evolution of an Arss cluster in an xfel pulse 
with hu> — 350 eV, / = lau, and T = 100 fs. Upper panel: 
Average time a quasi-free electron stays at a particular cluster 
atom. Middle panel: Energy range of quasi-free electrons 
(mean value ± standard deviation) compared to the average of 
the barriers between ions in the cluster. Dashed line: Radius 
of the cluster, see right axis. Lower panel: Average of ratios 
-Bfrco/Bbound and powers of it, see text. 
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FIG. 4: Average charge per atom for two cluster sizes Ari3 
and Arcs produced by an xfel pulse (same parameters as 
in Fig. |lj) as a function of the field strength /. Dot-dashed 
line: Restricted cluster calculation where intra-cluster screen- 
ing was precluded. 

ization due to the cluster space charge we have repeated 
our cluster calculations without the possibility for tun- 
neling of electrons between cluster ions. Figure ^ shows 
the final charges from these restricted calculations for 



Ari3 and Ar^ along with the full calculation for these 
clusters and the atom. The difference between the two 
cluster calculations with and without tunneling (marked 
by grey shading in Fig. ||) accounts for the derealization 
effect. The difference between the restricted cluster cal- 
culation (the dash-dotted line in Fig. ^) and that for the 
atom reveals the space charge effect. For field strengths 
/ > 0.3 au, where differences between atom and clusters 
appear, the space charge effect is initially weaker. This 
changes for stronger fields: Whereas for the smaller clus- 
ter Ari3 both are of the same magnitude at / = 10 au, for 
the larger cluster Ars5 the space charge effect dominates 
at this field strength. 

In conclusion, we have found suppression of ionization 
for small argon cluster compared to isolated atoms in the 
same pulse of intense X-ray radiation. This behavior is 
in striking contrast to that of rare-gas clusters in intense 
optical lasers. Two effects are responsible for the differ- 
ence: Firstly, the high positive space charge of the cluster 
hinders electron emission since the space charge is not 
compensated by a large quiver motion. Secondly, der- 
ealization of electrons in the cluster reduces photoioniza- 
tion as well as autoionization drastically. Both effects are 
more important for the larger cluster investigated and the 
relative weight of both effects depends on field strength 
and cluster size. Our findings indicate that in general the 
coupling of energy from the laser light to matter is less 
effective at high frequencies. This has important conse- 
quences for xfel imaging applications since it implies a 
higher damage threshold. That is, fragile samples might 
survive irradiation by intense pulses of high frequency 
better than anticipated. 

We would like to thank Christian Siedschlag for helpful 
discussions. 
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